Quantum Phase Transitions and the v = 5/2 Fractional Hall State in Wide Quantum Wells 
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We study the nature of the v = 5/2 quantum Hall state in wide quantum wells under the mixing of electronic 
subbands and Landau levels. We introduce a general method to analyze the Moore-Read Pfaffian state and its 
particle-hole conjugate, the anti-Pfaffian, under periodic boundary conditions in a "quartered" Brillouin zone 
scheme containing both even and odd numbers of electrons. We examine the rotational quantum numbers on 
the torus, and show spontaneous breaking of the particle-hole symmetry can be observed in finite-size systems. 
In the presence of electronic- subband and Landau-level mixing the particle-hole symmetry is broken in such a 
way that the anti-Pfaffian is unambiguously favored, and becomes more robust in the vicinity of a transition to 
the compressible phase, in agreement with recent experiments. 
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The quantized Hall state at the v = 5/2 Landau level 
(LL) filling factor [JJ has been the subject of significant re- 
cent interest due to a strong suspicion, with considerable sup- 
port from the numerical calculations J2-01, that it is described 
by the Moore-Read "Pfaffian" (Pf) state fij. This incom- 
pressible quantum fluid, which is a prototype state for non- 
Abelian exchange statistics fllzl], is found in the vicinity of 
the phase boundary with compressible phases characterized 
by stripe order and Fermi-liquid-like behavior yfl. Generally, 
the incompressible fluids of the fractional quantized Hall ef- 
fect (FQHE) ff8l— TToh possess protected gapless edge modes, 
and gapped bulk excitations which carry fractional charge and 
obey fractional statistics These attributes - quantiza- 

tion, fractionalization, and protection - represent the hallmark 
of topological phases Jl2ll . In the case of non-Abelian states, 
the degeneracy of the quasi-particle states may be suitable to 
implement a "fault- tolerant" quantum computation II 1 311 . 

The Moore-Read state, though defined in a half-filled LL, 
is not invariant under the particle-hole (P-H) transformation, 
and therefore its P-H conjugate partner - the "anti-Pfaffian" 
(APf) [14] - emerged as a competing candidate to describe 
the ground state of v = 5/2. In experiment, either Pf or 
APf is realized, depending on the explicit form of the P-H 
symmetry-breaking fields (e.g., LL mixing EMI). In the 
absence of those, the true ground state is selected by spon- 
taneous P-H symmetry breaking. A similar outcome should 
be reproducible in finite- size calculations, which have been 
known to capture remarkably well the fundamental aspects of 
FQHE physics ifioll . However, for technical reasons (see be- 
low), this never occurs for an even number of electrons which 
has been the assumption of most studies to date Jl8[[l9 1. 



Although Pf and APf have identical non-Abelian braiding 
properties in the bulk, they represent distinct phases of matter 
as reflected e.g., in their ed ge p hysics — a signature of the 
underlying topological order|jl2]. A number of recent exper- 
iments have focused on measuring the quasiparticle charge at 
v = 5/2 1E0U22I1 . and on detecting the non-Abelian statistics 
using edge-tunneling interferometry |23|,|24]]. These probes, 
while not definitive in all regards, are consistent with the non- 



Abelian statistics, and have provided additional insights into 
the nature of the ground state. In particular, the discovery 
of the counter-propagating mode [25J is consistent only with 
APf. 

At the same time, other experiments have focused on prob- 
ing the stability of v = 5/2 by driving the transitions from 
the incompressible to compressible phases, such as the Fermi 
liquid-like state igrl E7I1 and the anisotropic, stripe and ne- 
matic phases ll3l l28ll. This was done by tilting the mag- 
netic field JH], and by tuning the density in wide quantum 
well (WQW) samples J29[ J30Q. In the latter case, it was re- 
cently noticed [30] that the quantized v = 5/2 state becomes 
stronger in the vicinity of a transition to the Fermi liquid-like 
phase. 

In this Letter we introduce a new method to study the 
physics of the Moore-Read state and its P-H conjugate. We 
consider a compact torus geometry |3lll32|] with a "quartered" 
many-body Brillouin zone for both even and odd number of 
electrons N e . Regardless of the parity of N e , for the Moore- 
Read 3 -body Hamiltonian we obtain a zero-energy and zero- 
momentum ground-state, with bosonic ("magneto-roton") and 
fermionic ("neutral fermion" 113 311 ) collective modes at fixed 
N e . This is in stark contrast e.g. with the spherical geom- 
etry ll34ll . where the zero-energy ground-state only exists for 
even N e , and the fermionic mode can only be obtained for odd 
N e ll33ll . The essential physical similarity of the even and odd 
N e cases on the torus enables us to restrict to the latter case 
when Pf and APf can be further classified by their invariance 
under discrete rotations in high symmetry Bravais lattices. For 
special odd values of N e we recover spontaneous P-H symme- 
try breaking in finite systems as Pf and APf acquire different 
angular momenta compatible with periodic boundary condi- 
tions (PBCs), and become orthogonal to each other in finite 
systems. This formalism is applied to a realistic model of a 
wide quantum well with two subbands, SI (symmetric sub- 
band with n = 1 LL form-factor) and AO (antisymmetric sub- 
band with n = LL form-factor), where P-H symmetry is 
broken explicitly by the mixing of electronic subbands/LLs as 
a result of tuning the density [30]. We identify APf as the one 
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describing the ground state in these circumstances, and show 
that its gap increases prior to the transition to the compressible 
phase, in agreement with experiments 

We consider TV e electrons in a fundamental domain Li x L2 
subject to magnetic field Bz = V x A. An operator that 
translates a single electron and commutes with the Hamilto- 
nian is the magnetic translation operator which obeys a non- 
commutative algebra, leading to the quantization of the flux 
TV$ threading the system, z • (Li x L2) = 2ir£%N$, where 
£b = \JhjeB is the magnetic length. In a many-body sys- 
tem, symmetry classification is achieved by the help of the 
emergent many-body translation operators 1 32], which can be 
factorized into a center of mass and a relative part. The action 
of the former produces a characteristic degeneracy equal to q, 
where N e = pN, TV$ = qN, and TV is canonically assumed to 
be the greatest common divisor of TV e and TV$ such that p and 
q are coprime. The eigenvalue of the relative translation op- 
erator is a many-body momentum k 13211 that fully classifies 
the spectrum in the Brillouin zone TV x TV, with the excep- 
tion of high symmetry points where discrete symmetries may 
produce additional degeneracies, as we explain below. 

For the outlined algebraic derivation it is not essential 
that p, q be coprime numbers. In fact, enforcing this con- 
dition might hide the important physical features of a FQH 
state. This occurs for the Moore-Read state which possesses 
a pairing structure revealed in its fundamental root pattern 
110011001100 .. ., that defines the clustering properties on 
genus-0 surfaces l35ll . The corresponding filling factor v = 
1/2 should be viewed as v = 2/4 because the root pattern ad- 
mits two particles in each four consecutive orbitals. To incor- 
porate this clustering condition, we need to map the wavevec- 
tors onto a "quartered" Brillouin zone TV x TV 13^1 . We can 
view the quartered Brillouin zone as a result of "folding" the 
original TV x TV zone, Fig.[TJa). For TV even, the zone corner 
and midpoints of the zone sides all map to k = point (red 
points in Fig. Uta)). These wavevectors also define the sec- 
tors of the Hilbert space where a three-fold degenerate Moore- 
Read groundstate is obtained [6] . Moore-Read parent Hamil- 
tonian also possesses a zero-energy ground state for an odd 
number of electrons, which corresponds to a single unpaired 
electron with k = 0]. In this case, the folding still "com- 
pactifies" the original Brillouin zone, but the k = sector 
remains the same as before the folding. Thus, in the quartered 
Brilloiun zone, the ground-states of the Moore-Read Hamilto- 
nian are invariably obtained in the k = sector of the Hilbert 
space, as they should be for an incompressible liquid, and the 
3 -fold degenerate states are allowed to mix. 

The virtue of a quartered Brillouin zone becomes obvious 
when the full energy spectrum is studied as a function of mo- 
mentum, Fig. Hlb). Using the conventional definition of the 
Brillouin zone, there is no obvious structure in the low-lying 
excitation spectrum of the Moore-Read 3 -body Hamiltonian. 
However, if the same spectrum is replotted in a quartered 
zone, it reveals a clear bosonic mode (the "magneto-roton") 
and a fermionic mode ("neutral fermion") l33ll . Because of 
the Brillouin zone folding, in the even case (TV e = 14) we 




k x [2rc/L x ] 0.5 1 1.5 2 2.5 



FIG. 1 . (Color online) (a) An example of a squared many-body Bril- 
louin zone (open circles), and the mapping to a "quartered" Brillouin 
zone (black circles) for even TV e . Red circles denote sectors that map 
to the zero momentum k = 0. (b) Energy spectrum of the Moore- 
Read 3 -body Hamiltonian in a quartered Brillouin zone for TV e = 14 
and TV e = 13 (inset) particles. 



obtain three copies of the bosonic mode and a single copy of 
the fermionic one. In the odd number case (TV e = 13), the 
multiplicities are interchanged D6I1 . Contrary to spherical ge- 
ometry where spectra for different particle numbers have to be 
superimposed on the same plot to resolve the two modes, on 
the torus both modes are obtained for a fixed system size TV e . 
Additional advantage of PBCs is that one can access a quasi- 
continuum of the momenta k, as opposed to a much smaller 
subset of angular momenta on the sphere. This is achieved 
by adiabatic variation of the shape of the unit cell in terms 
of its aspect ratio | Li | / 1 L2 1 or the angle between vectors Li 
and L2, subject to a constraint that the area |Li x L2 1 remains 
fixed and equal to 2irl 2 B N<s>. In Fig. [T] we set the ratio equal to 
unity, and vary the angle between the square and the hexagon. 

Using a quartered Brillouin zone, we can also directly ad- 
dress the phase transition between the Moore-Read state and 
the composite Fermi liquid (CFL). Spherical geometry is not 
adequate for this purpose because the two states have differ- 
ent "shifts" l37ll . To capture the transition, we study the en- 
ergy spectrum of a Hamiltonian that interpolates between the 
model 3 -body Hamiltonian and the n = LL Coulomb in- 
teraction, A#3b + (1 — A)iJc- in a quartered Brillouin zone, 
the spectra for even and odd TV e are strikingly similar, and 
data in Fig. [2 corresponds to TV e = 11 which is also stud- 
ied in a different model below. The neutral fermion and the 
magneto-roton modes become significantly distorted and dif- 
ficult to identify for A < 1, nevertheless one can track their 
evolution until the eventual collapse of the gap for A = 0. 
At this point the ground state moves from k = to some 
k ~ ig 1 , and the system undergoes a second-order transition 
to the compressible phase. In contrast, for n = 1 LL Coulomb 
interaction, similar calculation does not lead to the gap clos- 
ing for any A. Adiabatic variation of PBCs also enables one to 
calculate e.g. the Hall viscosity, which is expected to diverge 
at A = point. Regime of small A might display an inter- 
esting crossover between type-I and type-II superconducting 
behavior iH. 

In order to address the competition between Pf and APf for 
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FIG. 2. (Color online) Transition between the Moore-Read and the 
CFL state. Energy spectrum of the Hamiltonian XHsh + (1 — X)Hc 
is plotted as a function of momentum for several A, illustrating the 
collapse of the neutral mode around k ~ i^ 1 . 



the ground state of generic (2-body) interactions, it is essential 
to consider all possible symmetries of the Hamiltonian at half 
filling. Particularly, the two symmetries of interest here are the 
P-H conjugation r p h, and discrete rotations compatible with 
PBCs ||39l|40|]. The first is represented by an anti-unitary op- 
erator whose effect is similar to the well-known case of time- 
reversal symmetry Eh. In the presence of both symmetries, 
extra (isolated) degeracies may occur which belong to con- 
jugate or co-representations of discrete rotations (32 



till]. In 



the case of even N e , rotational symmetry does not lead to any 
additional quantum numbers for the ground state at v = 5/2. 



For some excited states with k 



0, 



' P h 



-1 and these, 



following Kramer's theorem, are all doubly degerate. On the 
other hand, for odd N P and k =0, r r 



ph 



1 and only isolated 

degeneracies are possible. The ground state in this case is ei- 
ther unique or a doublet. For any geometry other than hexag- 
onal, the ground state is a singlet. For hexagonal geometry, 
the ground state is also a singlet if the number of electrons is 
given by iV e = 6 m + 1, where m G Z. For other N e , the 
ground state is a doublet B42I1 . 

To understand these trends for generic Hamiltonians it is 
helpful to consider the rotational properties of Pf and APf 
model states. On the torus with n-fold point symmetry, the 
angular momentum of APf, measured relative to Pf, is given 
by AM = 2N pa ir (mod n), where N pa i r is the number of 
paired electrons. For n = 2 or n = 4 (square or lower sym- 
metry), AM = since N pair is always even. Thus there is 
no symmetry reason for Pf and APf (being the eigenstates of 
different Hamiltonians) to be orthogonal. On the other hand, 
for n = 6, AM =^ if N e ^ 6m + 1, and the two states 
are necessarily orthogonal. It is under precisely these condi- 
tions that the doublet ground states are observed. These seem- 
ingly unrelated events are another confirmation that in the case 
of the Coulomb interactions the system is in the Moore-Read 
phase. There are always doublets present in the spectrum, but 



they describe the ground state only in cases where Pf and APf 
have AM / 0. Evidently, the rotational quantum numbers 
of the doublet match those of Pf and APf. There will be no 
spontaneous breaking of the P-H symmetry in the absence of 
such degeneracies until the thermodynamic limit is reached. 
Achieving this property for finite sizes makes the comparison 
of the exact ground state with Pf or APf much cleaner than the 
case of even N e , for example. As we show below, LL mixing 
splits the doublets in such a way that each member has a finite 
overlap with either Pf or APf, while having zero overlap with 
the other. 

We illustrate the ideas above on a model of the WQW in 
which electrons can populate two "active" subbands with LL 
indices n = and n = 1 114311 . We refer to these levels as AO 
and SI, where S,A stands for the wavefunction in the perpen- 
dicular z G [0,w] direction, w being the width of the well. For 
simplicity, we assume reflection symmetry around z = w/2, 
and the two subbands are given by symmetric/antisymmetric 
infinite square well wavefunctions ips = -\/2/wsm(7rz/w), 
(fA = sm(2ir z/w). The remaining subbands are either 

completely filled or completely empty, and excitations to them 
are forbidden. The energy splitting between the subbands, 
As as, can be realistically tuned by electrostatic gates 
but in our calculation it is assumed to be an independent pa- 
rameter. This model is expected to provide a realistic descrip- 
tion of a number of recent experiments on the WQWs where 
FQH states were probed by tuning the effective interaction via 
subband/LL mixing. In particular, we focus on the v = 5/2 
considering a half-filled S 1 level mixing with an AO level. It 
is assumed that the electron spin is fully polarized 

In Fig. Oa) we plot the neutral gap and mean value of the 



pseudospin operator S z = ^X^( c s,; c s,i _ c J A4 c A ,i) in the 
ground state for N e = 10 particles. We set w/£b = 2.7, 
which roughly agrees with the experimental width |28l [3(1 • 
When As as is large, excitations to AO level are costly, and 
the ground state is fully polarized in S 1 level and it is of the 
same nature as the one that is usually observed in wide sam- 
ples. As As as becomes smaller, the difference of n = 
and n = 1 LL form factors makes it increasibly favorable to 
promote particles into AO level and reduce the correlation en- 
ergy. Eventually, all particles migrate to AO subband, where 
they form a composite Fermi liquid; as shown in Fig. Oa), 
this happens slightly before the actual coincidence of the two 
subbands. The step-like behavior of (S z ) suggests the transi- 
tion to be very sharp, and we find it to be more affected by 
the difference in LL (rather than the subband) form-factors. 
Right before the transition, excitations from SI to AO lead to 
an increase of the neutral gap of the system, as measured ex- 
perimentally We believe the increase of the gap to be 
an intrinsic feature of the system, but limitations on the sys- 
tem sizes attainable by exact diagonalization prevent us from 
performing proper finite- size scaling of the gap. 

Finally, we investigate the nature of the ground state be- 
fore the transition to the compressible phase. In Fig. Ob) we 
compare the N e = 10 electron ground state (projected to SI 
level) with the Pf and APf wavefunctions defined on a hexag- 
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FIG. 3. (Color online) (a) Neutral gap and mean value of the pseudospin S z for iV e = 10 particles at half filling in a WQW (w/Ib = 2.7). 
Transition is driven by changing As as, and relative position of the electronic subbands is schematically shown in the inset, (b) Overlap 
between the exact iV e = 10 ground state (projected onto SI subband) and the Pf/APf wavef unctions. In the transition region (As as ~ 0.028) 
APf becomes favored, (c) Energy spectrum for infinite As as (black diamonds) with an exact doublet at k = 0, and the splitting of the doublet 
when Asas is reduced (red triangles), (d) Same as in (b) but for N e = 11, when Pf and APf are mutually orthogonal. Asas breaks PH 
symmetry and unambiguously selects APf (95% overlap) over Pf (zero overlap). Inset shows the effect of varying the V\ pseudopotential. 



onal unit cell. At large Asas, the PH symmetry is preserved, 
and Pf/APf have identical overlaps with the exact state. In the 
transition region, PH symmetry is lifted, and APf overlap is 
growing while that of the Pf is decreasing. However, Pf and 
APf have a signicant overlap with each other for N e = 10 
particles 11711 . To avoid this problem, we instead consider 
N e — 1 1 case, when discrete symmetry forces Pf and APf to 
be orthogonal. At the same time, the Coulomb ground state, 
which is exactly twofold degenerate for infinite Asas, splits 
when Asas is reduced (Fig.Oc)). The size of the full Hilbert 
space for N e = 11 is too large to be directly diagonalized, 
so we follow the approach of limiting the number of excita- 
tions to the higher subband H 1 Til . The convergence is found to 
be rapid, and the essential properties of the ground state and 
the transition point can be captured accurately by allowing for 
not more than 3 electrons in a higher subband. The effect of 
a small amount of symmetry-breaking by Asas is sufficient 
to select APf with 95% overlap with the ground state, while 
the Pf overlap drops to zero (Fig.^d)). The APf overlap can 
be increased further by varying the V\ pseudopotential in a 
"SU(2)-invariant" manner II 1 Tfl (Fig. Efd), inset). These re- 
sults provide unambiguous evidence that the ground state of 
the WQW is described by APf. 

In summary, we have demonstrated that PBCs for special 
finite systems with an appropriately defined Brillouin zone 
are essential for the study of the non-Abelian FQH states, 
their collective excitation spectra, and the effects of particle- 
hole symmetry breaking in WQWs under LL/subband mix- 
ing. Although we have presented results for the case of SI /AO 
mixing, the more usual case of S1/S2 mixing (applicable to 
narrower quantum well samples) selects APf to describe the 
ground-state via the same mechanism. This conclusion is 
not affected by the mixing of even higher LLs (S3 and be- 
yond). The sharpness of the transition and the "population 
inversion" illustrated in Fig.[3]might be relevant for the emer- 
gence of lowest LL physics under extreme tilting of the mag- 
netic field Q. 
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